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Abstract: Traditionally, longitudinal shower profiles are reconstructed in fluorescence light 
experiments by treating the Cherenkov light contribution as background. Here we will argue 
that, due to universality of the energy spectra of electrons and positrons, both fluorescence 
and Cherenkov light can be used simultaneously as signal to infer the longitudinal shower 
development. We present a new profile reconstruction method that is based on the analytic 
least-square solution for the estimation of the shower profile from the observed light signal 
and discuss the extrapolation of the profile with a Gaisser-Hillas function. 
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Introduction 

During its passage through the atmosphere 
of the earth an extensive air shower excites 
nitrogen molecules of the air, which subse- 
quently radiate isotropically ultraviolet fluo- 
rescence light. Since the amount of emitted 
light is proportional to the energy deposited, 
the longitudinal shower development can be 
observed by appropriate optical detectors such 
as HiRes [1], Auger [2] or TA [3]. 
As part of the charged shower particles travel 
faster than the speed of light in air, Cherenkov 
light is emitted in addition. Therefore, in gen- 
eral a mixture of the two light sources reaches 
the aperture of the detector. 
In the traditional method [4] for the recon- 
struction of the longitudinal shower develop- 
ment the Cherenkov light is iteratively sub- 
tracted from the measured total light. The 
drawbacks of this ansatz are the lack of con- 
vergence for events with a large amount of 
Cherenkov light and the difficulty of propagat- 
ing the uncertainty of the subtracted signal to 
the reconstructed shower profile. 
It has already been noted in [5] that, due to 
the universality of the energy spectra of the 
secondary electrons and positrons within an 
air shower, there exists a non-iterative solution 



for the reconstruction of a longitudinal shower 
profile from light detected by fluorescence tele- 
scopes. 

Here we will present the analytic least-square 
solution for the estimation of the shower profile 
from the observed light signal in which both, 
fluorescence and Cherenkov light, are treated 
as signal. 

Scattered and Direct Light 

The non-scattered, i.e. direct fluorescence light 
emitted at a certain slant depth X$ is measured 
at the detector at a time tj . Given the fluores- 
cence yield Y? [6, 7] at this point of the atmo- 
sphere, the number of photons produced at the 
shower in a slant depth interval AXi is 

N*(Xi) — Yl Wi AXi, 

where Wi denotes the energy deposited at slant 
depth Xi (cf. Fig. 1). These photons are 
distributed over a sphere with surface 4 ir rf , 
where r, denotes the distance of the detector. 
Due to atmospheric attenuation only a frac- 
tion Ti of them can be detected. Given a light 
detection efficiency of e, the measured fluores- 
cence light flux y\ can be written as 



V, 



diYfwi AXi, 



(1) 
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Figure 1: Illustration of the isotropic fluorescence light emission (circles), Cherenkov beam along 
the shower axis and the direct (left) and scattered (right) Cherenkov light contributions. 



where the abbreviation d; = f Ti -> was used. 
For the sake of clarity the wave length depen- 
dence of Y, T and e will be disregarded in the 
following but be discussed later. 
The number of Cherenkov photons emitted at 
the shower is proportional to the number of 
charged particles above the Cherenkov thresh- 
old energy. Since the electromagnetic compo- 
nent dominates the shower development, the 
emitted Cherenkov light, N^, can e calculated 
from 

N°(X i )=YpN?AX i , 

where N° denotes the number of electrons and 
positrons above a certain energy cutoff, which 
is constant over the full shower track and not to 
be confused with the Cherenkov emission en- 
ergy threshold. Details of the Cherenkov light 
production like these thresholds are included 
in the Cherenkov yield factor [5, 8, 9, 10]. 

Although the Cherenkov photons are emitted 
in a narrow cone along the particle direction, 
they cover a considerable angular range with 
respect to the shower axis, because the charged 
particles are deflected from the primary parti- 
cle direction due to multiple scattering. Given 
the fraction fc(Pi) of Cherenkov photons emit- 
ted at an angle Pi with respect to the shower 
axis [8, 10], the light flux at the detector aper- 
ture originating from direct Cherenkov light is 

y? d = d l f c (p i )Y z c AX z N^. (2) 

Due to the forward peaked nature of 
Cherenkov light production, an intense 



Cherenkov light beam can build up along 
the shower as it traverses the atmosphere 
(cf. Fig. 1). If a fraction / s (/3,) of the beam 
is scattered towards the detector it can con- 
tribute significantly to the total light received. 
In a simple one-dimensional model the number 
of photons in the beam at depth Xi is just 
the sum of Cherenkov light produced at all 
previous depths Xj attenuated on the way 
from Xj to Xi by Tji'. 

i 

N^ m (Xi) = Tji Y f ax j N j- 

3=0 

Similar to the direct contributions, the scat- 
tered Cherenkov light received at the detector 
is then 

i 

y? 8 = ./, fM J2 7 ," Y f AX 3 (3) 

3=0 

Finally, the total light received at the detec- 
tor at the time ti is obtained by adding the 
scattered and direct light contributions. 

Shower Profile Reconstruction 

The aim of the profile reconstruction is to es- 
timate the energy deposit and/or electron pro- 
file from the light flux observed at the detector. 
At first glance this seems to be hopeless, since 
at each depth there are the two unknown vari- 
ables Wi and N? , and only one measured quan- 
tity, namely yi. Since the total energy deposit 
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is just the sum of the energy loss of electrons, 
Wi and Nf are related via 



Wi 



N? [°° f e (E,Xi)w e (E)dE, 
Jo 



(4) 



where f e (E,Xi) denotes the normalized elec- 
tron energy distribution and w e (E,Xi) is the 
energy loss of a single electron with energy E. 
As it is shown in [9, 5, 10], the electron energy 
spectrum f c (E,Xi) is universal in shower age 
Si — 3/(1 + 2X max /Xi), i.e. it does not depend 
on the primary mass or energy, but only on 
the relative distance to the shower maximum, 
I max . Eq. (4) can thus be simplified to 

Wi = Nf on. 

where on is the average energy deposit per elec- 
tron at shower age Sj. With this one-to-one 
relation between the energy deposit and the 
number of electrons, the shower profile is read- 
ily calculable from the equations given in the 
last section. For the solution of the problem, 
it is convenient to rewrite the relation between 
energy deposit and light at the detector in ma- 
trix notation: Let y = (j/i, 1/2, ■ ■ ■ , y n ) T be the 
n-componcnt vector (histogram) of the mea- 
sured photon flux at the aperture and w = 
(u>i, W2, ■ ■ ■ , w n ) T the energy deposit vector at 
the shower track. Using the ansatz 



C • w 



(5) 



the elements of the Cherenkov-fluorescence 
matrix C can be found by a comparison with 
the coefficients in equations (1), (2) and (3): 



Cij — < 



0, 

4 + 



c ■ 



i < j 



i > 3, 



(6) 



where 



4 = di{Y? + fc(PJYf/ai) AX, 



and 



cl^d,m)T ri Y 3 c /a,AX r 

The solution of Eq. (5) can be obtained by in- 
version, leading to the energy deposit estima- 
tor w: 

W = C _1 -y. 



Due to the triangular structure of the 
Cherenkov-fluorescence matrix the inverse can 
be calculated fast even for matrices with large 
dimension. As the matrix elements in (6) are 
always > 0, C is never singular. 
The statistical uncertainties of w are obtained 
by error propagation: 

v w = c- 1 v y (c T ) _1 . 

It is interesting to note that even if the 
measurements are uncorrelated, i.e. their 
covariance matrix V y is diagonal, the calcu- 
lated energy loss values Wi are not. This is, 
because the light observed during time interval 
i does not solely originate from u>i, but also 
receives a contribution from earlier shower 
parts Wj, j < i, via the 'Cherenkov beam'. 



Wavelength Dependence 

Until now it has been assumed that the shower 
induces light emission at a single wavelength A. 
In reality, the fluorescence yield shows distinct 
emission peaks and the number of Cherenkov 
photons is proportional to p-. In that case, 
also the wavelength dependence of the detec- 
tor efficiency and the light transmission need 
to be taken into account. Assuming that a 
binned wavelength distribution of the yields is 
available (Y,fe = J_^ fe J^^ 1^(A) dA), the above 
considerations still hold when replacing cf and 
4 in Eq. (6) by 



cf = AXi ]T d lk {Yi + / c (ft) Yg/ai) 



and 



c? j =AX j '£ i d ik f a (l3 i )T Jik Y%/a j 



where 



The detector efficiency and transmission co- 
efficients Tik and 7^ are evaluated at the 
wavelength Afe. 
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Shower Age Dependence 

Due to the age dependence of the electron 
spectra f e (E,Si), the Cherenkov yield factors 
Yp and the average electron energy deposits 
on depend on the shower maximum, which is 
not known before the profile has been recon- 
structed. Fortunately, these dependencies are 
small: In the age range of importance for the 
shower profile reconstruction (s e [0.8, 1.2]) a 
varies only within a few percent [10] and Y c by 
less than 15% [5]. Therefore, a good estimate 
of a and Y c can be obtained by setting s = 1. 
After the shower profile has been calculated 
with these estimates, A max can be determined 
and the profiles can be re-calculated with an 
updated Chercnkov-fluorescence matrix. 



their average values (X ) and (A). The new 
minimization function is then the modified \ 2 



2 2 , (-^0 - (Af)}) 

X =Xgh + 



V; 



+ 



(A -(A)) 2 



where the variance of Xq and A around their 
mean values are in the denominators. 
In this way, even if Xgh i s n °t sensitive to Xq 
and A, the minimization will still converge. 
On the other hand, if the measurements have 
small statistical uncertainties and/or cover 
a wide range in depth, the minimization 
function is flexible enough to allow for shape 
parameters differing from their mean values. 
These mean values can be determined from 
air shower simulations or, preferably, from 
high quality data profiles which can be recon- 
structed without constraints. 



Gaisser-Hillas Fit 

The knowledge of the complete profile is re- 
quired for the calculation of the Cherenkov 
beam and the shower energy. If due to the 
limited field of view of the detector only a part 
of the profile is observed, an appropriate func- 
tion for the extrapolation to unobserved depths 
is needed. A possible choice is the Gaisser- 
Hillas function [11] which was found to give 
a good description of measured longitudinal 
profiles [12]. It has only four free parameters: 
A max , the depth where the shower reaches its 
maximum energy deposit w max and two shape 
parameters Xq and A. 

The best set of Gaisser-Hillas parameters p can 
be obtained by minimizing the error weighted 
squared difference between the vector of func- 
tion values fcH and x, which is 

X 2 GH = [w-f(p)] T Vw 1 [w-f(p)] 

This minimization works well if a large frac- 
tion of the shower has been observed below and 
above the shower maximum. If this is not the 
case, or even worse, if the shower maximum is 
outside the field of view, the problem is under- 
determined, i.e. the experimental information 
is not sufficient to reconstruct all four Gaisser- 
Hillas parameters. This complication can be 
overcome by weakly constraining Xq and A to 
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